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Analysis of Vortex Sheets 

Sijue Wu* 
Abstract 

We consider the motion of the interface separating two domains of the 
same fluid that moves with different velocity along the tangential direction of 
the interface. We assume that the fluids occupying the two domains are of 
constant densities that are equal, are inviscid, incompressible and irrotational, 
and that the surface tension is zero. We discuss results on the existence and 
uniqueness of solutions for given data, the regularity of solutions, singularity 
formation and the nature of solutions after the singularity formation time. 
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1. Introduction 

Vortex dynamics is of fundamental importance for a wide variety of concrete 
physical problems, such as lift of airfoils, mixing of fluids, separation of boundary 
layers, and generation of sounds. In mathematical analysis, one often neglects 
surface tension and viscosity, when they are small in the real physical problem. 
This necessitates justifying such simplifications. 

In this paper, we consider the motion of the interface separating two domains 
of the same fluid in R 2 that moves with different velocity along the tangential 
direction of the interface. We assume that the fluids occupying the two domains 
separated by the interface are of constant densities that are equal, are inviscid, 
incompressible and irrotational. We also assume that the surface tension is zero, 
and there is no external forces. The interface in the aforementioned fluid motion is 
a so called vortex sheet. We want to study the following problem: 

Given a vortex sheet initial data, is there a unique solution to this problem? 
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In general, there are two approaches to the aforementioned problem. One is 
to solve the initial value problem of the incompressible Euler equation in R 2 : 



where the initial incompressible velocity v G Lf oc (R 2 ), in which the vorticity uj — 
cvltIv is a finite Radon measure. Here v is the fluid velocity p is the pressure, 
and the density of the fluid is assumed to be one. Notice that a vortex sheet 
gives a measure valued vorticity supported on the interface. This approach was 
posed by DiPerna and Majda in 1987 [9]. In 1991, J.M. Delort [8] proved the 
existence of weak solutions global in time of the 2-D incompressible Euler equation 
(1) for measure- valued initial vorticity in H^ ].{R 2 ) that has a distinguished sign. 
However the problem of uniqueness of the weak solution is still unresolved. In 
1963, Yudovich [33] obtained the existence and uniqueness of weak solutions of the 
2-D incompressible Euler equation (1) for initially bounded vorticity. The best 
results on uniqueness upto date are given by Yudovich [34] and Vishik [31] for weak 
solutions with vorticity in a class slightly larger than L°° . This does not include 
vortex sheets, which admit measure-valued vorticity. Examples of weak solutions 
with the velocity field v G L 2 (R 2 x (— T, T)) that is compactly supported in space- 
time was constructed by V. Scheffer [29] and later by A. Shnirelman [30]. This gives 
non-uniqueness of weak solutions in L 2 (R 2 x (— T, T)). However non-uniqueness in 
the physically relevant class of conserved energy v G L°°([0, oo), Lf oc (R 2 )) remains 
open. Numerical evidences of non-uniqueness of weak solutions for vortex sheet 
data can be found in [25], [18]. 

Furthermore, weak solutions give little information of the specific nature of the 
vortex sheet evolution. For instance, does the vorticity remain supported on a curve 
for a later time given that the initial vorticity is supported on a curve in R 2 ? Assume 
further that the free interface between the two fluid domains remains a curve in R 2 
at a later time, equation (1) can be reduced to an evolutionary differential-integral 
equation along the interface. This is the Birkhoff-Rott equation, written explicitly 
by Birkhoff in [2] and implied in the work of Rott [26]. The second approach uses 
the Birkhoff-Rott equation as a model for the evolution of the vortex sheet. 



For convenience, we use complex variable z = x + i y to denote a point in R 2 . 
~z = x — iy denotes the complex conjugate and f x = d x f is the partial derivative of 
the function /. H s indicates Sobolev spaces. 

In search of the equation for the evolution of the vortex sheet, we suppose that 
at time t > the vorticity is a measure supported on the curve T(t) given by the 
complex position £ = £(s, t) in the arclength s, in which £(0, t) is the particle path 
of a reference particle; and on this curve the vorticity density is 7 = 7(s, t). That 





(1) 



2. The Birkhoff-Rott equation 
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is the vorticity at time t is uj(x, y, t) satisfying 

JJ <f>(x,y)u)(x,y,t)dxdy = J ^{s,t))j{s,t)ds, for any <f> e C^(R 2 ). 
From the Biot-Savart law, the velocity field v induced by the vorticity is given by 

Notice that the velocity is discontinuous just on T(t). We define the velocity on the 
sheet as the average of the velocities at the two sides of the sheet, that is given by 
the principle value integral: 

*«'•'>•'> = dW {(J) ( -'tt.o '*- (2) 

As suggested by the properties of the Euler equation, we assume that the vortex 
sheet is convected by the average velocity (2), and the vorticity is conserved along 
the particle path. We arrive at the evolution equation of the vortex sheet: 

t) + a(s, t) = v{£{8, t),t) 

lt {s,t)+d s (a(s,t)j( S ,t)) = 

where a(s, t) is a real valued function satisfying a(0, t) = 0. 

A rigorous justification of the equivalence between equation (3) and equation 
(1) for smooth graphs £(s,i) and smooth vortex strength j(s,t) can be found in 
[19]. It is not hard to extend it to all smooth curves. 

Assume a(s, t) = J Q S 7(s', t) ds' defines an increasing function of s, we make a 
change of variables: z(a, t) = £(s(a,t),t), in which s(a,t) is the inverse of a(s,t): 
a(s(a,t),t) = a. We get from equation (3) the Birkhoff-Rott equation 

d t z(a, t) = -L p. v . f 1 dp. (4) 

2ni J z(a, t) - z(p, t) 

Notice that z = z(a,t) is a parameterization of the vortex sheet in the circulation 
variable a, l/|z Q | = 7 is the vortex strength. A steady solution of (4) is the flat 
sheet z — a. 

Equation (4) has been under active investigations over the last four decades. 
A well-known property of (4) is that perturbations of the flat sheet grow due to the 
Kelvin-Helmholtz instability, following from a linearization of equation (4) about the 
flat sheet. For given analytic data, Sulem, Sulem, Bardos and Frisch [28] established 
the short time existence and uniqueness of solutions in analytic class for 2-D and 
3-D vortex sheet evolution. Duchon and Robert [10] obtained the global existence 
of solutions of equation (4) for a special class of initial data that is close to the flat 
sheet. However, numerous results show that a vortex sheet can develop a curvature 
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singularity in finite time from analytic data. D.W. Moore [21] was the first to 
provide analytical evidence that predicts the occurrence and time of singularity 
formation, which was verified numerically by Meiron, Baker and Orszag [20] and 
by Krasny [13]. Caflisch and Orellana [3] proved existence almost up to the time 
of expected singularity formation for analytic data that is close to the flat sheet. 
Duchon and Robert [10] and Caflisch and Orellana [4] constructed specific examples 
of solutions of equation (4) where a curvature singularity develops in finite time 
from analytic data. The example of Caflisch and Orellana [4] has the form z(a, t) — 
a + S(a, t) + r(a, t), where 

S{a,t) = e(l - i){(l - e -*/2-i«)i+P _ (! _ e -t/2+ la y+^ 

is a solution of the linearized equation in which e is small, fi > 0; r(a, t) is the 
correction term that is negligible relative to S(a, t) in the sense that S(a, t) +r(a, t) 
exhibits the same kind of behavior as S(a, t) [4]. Notice that S(a,t) is an analytic 
function for t > 0, but S(a,0) has an infinite second derivative at a = for 
H E (0, 1). In fact, the (1 + i^)th derivative of S(a,0) for v > [i becomes infinite at 
a = 0. Now inverting time gives an example z(a, t) that is analytic at t < 0, but 
has an infinite second derivative at a = 0, t = 0. At the singularity formation time 
t = 0, the vortex strength l/\z a \ of this example satisfies 

< c < l/\z a \ < C < oo (5) 

for some constants c and C; and z(a, t) G C 1+P (R x [i , 0]) for < p < /i. 

These examples also show that the initial value problem of the Birkhoff-Rott 
equation (4) is ill-posed in C 1+l/ (R), v > 0, and in Sobolev spaces H S (R), s > 3/2 
in the Hadamard sense [4], [10]. Ill-posedness was also proved by Ebin [11] using a 
different approach. However the existence of solutions in spaces less regular than 
C 1+U {R) or H S (R), and the nature of the vortex sheet at and beyond the singularity 
time remained unknown analytically in general. 

This suggests that we look for solutions of the Birkhoff-Rott equation in the 
largest possible spaces where the equation makes sense. For the purpose of this 
paper, we consider functions z(a, t) so that for each fixed time t, both sides of the 
equation (4) are functions locally in L 2 , and on which the L 2 -analysis is available. 
This leads us to consider chord-arc curves, thanks to the work of G. David [7]. 

Another reason that chord-arc curves are to be considered is due to the nu- 
merical calculation of Krasny [14] [15]. Krasny studied the evolution of the vortex 
sheet beyond singularity using the vortex blob method. He found that the approx- 
imating solutions have the form of a spiral beyond singularity. Convergence of the 
approximating sequence to a weak solution of the Euler equation (1) was proved 
by J-G. Liu and Z-P. Xin [17], under the assumption that the initial vorticity has a 
distinguished sign. A special example of chord-arc curves is a logarithmic spiral. 
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3. Chord-arc curves and some recent results 

Let r be a rectifiable Jordan curve in R 2 given by £ = £(s) in the arclength s. 
We say T is a chord-arc curve, if there is a constant M > 1, such that 



The infimum of all such constants M is called the chord-arc constant. 

For a chord-arc curve £ = £(s), s the arclength, it is proved in [6] that £'(s) 
exists almost everywhere, and there is a choice of the argument function b € BMO, 
with £'(s) = e lb ^ . In particular, if the chord-arc constant is close to 1, there is 
a choice of b € BMO, such that ||6|| bmo is close to 0. Moreover the subset of all 
those functions b is an open subset of BMO. And if & £ BMO, and ||6||smo < 1> 
£(s) = £o + Jo e J ^ s ) ds' defines a chord-arc curve. 

Examples of chord-arc curves include Lipschitz curves and logarithmic spirals 
r = ±e 9 , 9 e R, where (r, 0) is the polar coordinates. 

A Theorem of G. David [7] states that 

Theorem (G. David [7]). For all chord-arc curves T : £ = £(s), s £/ie arclength, 
the corresponding Cauchy integral operator Cr, where 



is bounded from L 2 (ds) to L 2 {ds). 

In fact, the result of G. David [7] is stronger than stated above. He proved that 
the Cauchy integral operator Cr is bounded from L 2 (ds) to L 2 (ds) if and only if T 
is a regular curve. A rectifiable curve T is said to be regular if there is a constant 
M such that for every r > and every disc D with radius r, the length of T n D 
does not exceed Mr. A chord-arc curve is regular but not vice versa. 

Now we go back to the Birkhoff-Rott equation. Notice that the Biot-Savart 
integral representing an incompressible velocity field v in terms of the vorticity u> 
may be divergent if u> does not vanish fast enough at infinity, even if the velocity 
field v is well defined, we extend the definition of the Birkhoff-Rott equation (4) by 
considering the differences of the velocities between any two points: 



1*1 -* 2 | <M|£( Sl )-£(s 2 )| 



for all s\, s 2 . 




, _ J_ f r 1 1 

[a, t) z t {a , t) - ^ p.v. J^J ^ t) t) - zK t) ^ t) 



}d(3 



(6) 




2?ri P ' ' J m>N {z{a, t) - z(J3, t))(z{a', t) - z{/3, t)) 



for all (a,t), (a',t), and some N > |a| + |a'| + 1. This admits a larger class 
of solutions. In particular, the integral on the right hand side of equation (6) is 
convergent for those similarity solutions considered in [12], [23]- [25], which otherwise 
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give divergent Cauchy integrals in (4) due to the divergent contributions from the 
vorticities at infinity. It follows from the Theorem of G. David that the integral on 
the right hand side of the equation (6) is convergent for a.e. (a,t), (a',t) and is 
in L°°([0, T], Lf oc (da) x Lf oc (da')) for the solutions considered in Theorem 1 in the 
following. 

Roughly speaking, if a function z = z(a,t) satisfies equation (4), it will also 
satisfy equation (6). On the other hand, if z = z(a, t) satisfies (6), and the Cauchy 
integral 



1 f 1 

liri P J z(a, t) — - 



2tu J z(a,t) - z(/3,t) 

is convergent, then there is a function c = c(t), such that z(a,t) + c(t) satisfies 
equation (4). 

For a local integrable function / = /(a) defined on (a, 6), we say / is of 
bounded local mean oscillation on (a, b) if there exists So > such that 

||/||bmo(o,6),« = SU P T7T / !/(") - fi\ da < °°> 

all IC(a,b),\I\<So I 1 I J I 

here // = |jy j I f{a)da, I is an interval. We say / is analytic on (a, b) if / 6 
C°°(a,b), and for any compact subset K of (a, 6), there is a constant p > 0, such 
that 



Efe / ISS7(a)| 2 da<oo. 



m—o 



Notice that lnz is multi- valued for complex number z. In the following, \nz a refers 
to one choice of the multi-values. We have the following results concerning the 
solutions of the Birkhoff-Rott equation (6) (or (4)). 

Theorem 1 [32]. Assume that z G i^QO, T], L? oc (R)) n L 2 ([0, T], Hl oc {R)) is a 
solution of the Birkhoff-Rott equation (6) for <t <T, satisfying that 

1. There are constants m > 0, M > 0, independent oft, such that 

m\a-f3\<\z(a,t)-z(f3,t)\<M\a-[3\ for all a, [3, < t < T. (7) 

2. For the interval (a,b), there exists So > 0, independent oft, such that for 
allO<t<T, 

\\lnz a {-,t)\\ BM o(a,b),So <C{m,M), (8) 

here C(m,M) is a universal constant depending on m and M. Assume further that 
lnz Q e L 2 ([0,T],Lf oc (R)). Then z a e C{(a,b) x (0,T)), and for each t G (0,T), 
^a(')io) is analytic on (a,b). 

Remark Notice that assumption 2. is satisfied if In z a G C([a, b] x [0, T]). Assump- 
tion 1. is equivalent to assuming that the vortex strength 7 = l/\z a \ is bounded 
away from and 00, with bounds independent of t, and z = z(-,t) defines a chord- 
arc curve for each fixed t G [0,T], with the chord-arc constant independent of t. 
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Assumption 1. can be relaxed by requiring that l/\z a (a, t)\ is bounded away from 
and oo for a G (a, b) only, and by assuming some weaker conditions for a G" (a, b). 

A version of Theorem 1 under assumptions that the solution of the Birkhoff- 
Rott equation (4): z = z(a,t) G C 1+Po , for some p > 0, the vortex strength 
Co < 7(a, i) < Co for some constants Co > 0, Co > 0, and the vortex sheets T(t) are 
closed Jordan curves is also obtained by Lebeau [16] using an independent approach. 

As a consequence of Theorem 1, the example constructed by Caflisch and 
Orcllana [4] will fail to satisfy properties 1 and 2 as stated in Theorem 1 near the 
singularity after the singularity formation time. 

Let Rez and Imz be the real and imaginary parts of the complex number z 
respectively. Regarding the existence of solutions, we have the following 

Theorem 2 [32]. For any real valued function wq G Hi{R), there exists T — 
T{ \\vjq\\ 3 ) > 0, such that the Birkhoff-Rott equation (6) has a solution z = z(a,t) 

for < t < T, satisfying \nz a G L°°([0, T], Hi (R)) n Lip([0, T], (R)) and 
Im{(l+i) lnz a (a, 0)} = wo(a), with the properties that there exist constants m > 0, 
M > independent oft, such that 

m\a -/3\< \z(a,t) - z{[3,t)\ < M\a - [3\, for all a, /3, < t < T; 

and there exists So > 0, independent of t, such that 

\\lnz a {;t)\\ BM o(R),6 <C(m,M), for < t < T, 

here C(m,M) is the universal constant as in Theorem 1. 

Theorem 2 states that if only half of the data z a (a, 0) is given, there is a 
solution of the Birkhoff-Rott equation (6) for a finite time period. Theorem 2 is a 
generalization of the existence result of Duchon and Robert [10] to general data. 

The following result implies that Theorem 2 is optimal, in the sense that in 
general, there is no solution of the Birkhoff-Rott equation satisfying properties 1. 
and 2. as stated in Theorem 1 beyond the initial time t = for arbitrarily given 
data. 

Theorem 3 [32]. Assume that z G /f 1 ([0, T], Lf oc (R)) n L 2 {[0, T], Hl oc {R)) is a 
solution of the Birkhoff-Rott equation (6) for 0<t<T,T>0, satisfying the 
property 1. and property 2. on some interval (a, b) as stated in Theorem 1. Assume 
further that lnz a G L 2 ([0,T], Lf oc (Rj), and Wo — Im{(l + z)lnz a (-,0)} is analytic 
on (a,b). Then z a G C((a,b) x [0,T)) and Re{(l + i) lnz Q (-, 0)} is also analytic on 
(a,b). ' 

4. Open questions 

The reason that Theorem 1-3 holds is that under the assumptions 1. and 2. 
in Theorem 1, the Birkhoff-Rott equation (6) is of "elliptic" type on (a, b) x [0, T]. 
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It would be interesting to sec whether the assumption 2 that requires small local 
mean oscillation on lnz Q (-, t) can be removed, or to see whether one can construct 
a non-smooth solution of equation (6) that violates (8). A good place to start is to 
construct similarity solutions of the Birkhoff-Rott equation. Similarity solutions are 
studied numerically in Pullcn [23] [25] , Pullen and Phillips [24] , and analytically in 
Kambc [12]. However, the similarity solutions in [12], [23]- [25] violate (7). It would 
also be interesting to relax the assumption 1. by considering vortex strength that 
is not necessarily bounded away from and infinity on the interval (a, b). 

Theorem 1-3 implies that a solution of the Birkhoff-Rott equation is either 
analytic, or doesn't satisfy properties 1. and 2. in Theorem 1. And there are in 
general no solutions of reasonable regularity (properties 1. and 2.) beyond the 
initial time for an arbitrarily given data. Are there solutions of the Birkhoff-Rott 
equation in spaces of even less regularity? Notice that we can define the solutions 
of the Birkhoff-Rott equation in the distribution sense by 

')*(«)**) =h IS ztt- Mi) ****' f ° r aU " E °~ 

This admits some even larger classes of solutions z = z(a,t). Is there a solution in 
the distribution sense for the Birkhoff-Rott equation for any given data? How is it 
related to the weak solution of the Eulcr equation? 

Theorem 1-3 might as well suggest that the vortex sheet in general fails to be 
a curve beyond the initial time for general data. Therefore it becomes interesting 
to study the vortex layers or considering the effects of viscosity. Numerical analysis 
of the vortex layers can be found in Baker and Shelley [1] etc. 
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